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Neutral pion transverse momentum (pt) spectra at midrapidity (\y\ < 0.35) were measured in 
Cu+Cu collisions at y'sjvjv = 22.4, 62.4, and 200 GeV. Relative to n° yields in p+p collisions scaled 
by the number of inelastic nucleon-nucleon collisions (JV co ii) the 7r° yields for px <^ 2 GeV/c in 
central Cu+Cu collisions are suppressed at 62.4 and 200 GeV whereas an enhancement is observed 
at 22.4 GeV. A comparison with a jet quenching model suggests that final state parton energy- 
loss dominates in central Cu+Cu collisions at 62.4 GeV and 200 GeV, while the enhancement at 
22.4 GeV is consistent with nuclear modifications in the initial state alone. 

PACS numbers: 25.75.Dw 



The measurement of particle yields at high transverse 
momentum (px ^ 2 GeV/c) has played a key role in char- 
acterizing the medium created in nucleus-nucleus colli- 
sions at the Relativistic Heavy Ion Collider (RHIC) [ILQ. 
Hadrons produced at sufficiently high px result from the 
interaction of quarks and gluons with high momentum 
transfer ("hard scattering") which can be described by 
perturbative quantum-chromodynamics (pQCD). These 
hadrons are produced as particle jets in the fragmenta- 
tion of the scattered partons. A scattered parton prop- 
agating through a quark-gluon plasma, a thermalized 
medium in which quarks and gluons are not confined 
in hadrons, loses energy ("jet-quenching") resulting in 
hadron yields at high px being suppressed [3j. Such a 
suppression was indeed observed in central Au+Au col- 
lisions at ytspfN = 130 and 200 GeV at RHIC, providing 
evidence for large color- charge densities in these systems 

BS0. 

Characteristic properties of the suppression of hadrons 
at high-px, e.g., the dependence on px and centrality, 
were studied in detail in Au+Au collisions at yj snn — 
200 GeV [5j. However, the energy dependence of hadron 
production in A+A collisions as predicted by jet quench- 
ing models @, H, is not well constrained by mea- 
surements. Work in this direction was presented in 
To study the energy dependence of jet- 



10, 11, LL 



quenching it is desirable to measure identified particles 
in the same colliding system over a large y/s^N range 
and to compare to p+p reference data measured in the 
same experimental setup. Identified particles provide an 
advantage over unidentified hadrons in that the inter- 
pretation is not complicated by the different contribu- 
tions from baryons and mesons. The study of Cu+Cu 



collisions is particularly useful because hadron suppres- 
sion in Au+Au collisions is observed for rather periph- 
eral collisions with a number of participating nucleons of 
A/part ~ 50 — 100 This N part range can be studied 
with reduced uncertainties in iV co ii with the smaller Cu 
nucleus. 

A critical parameter in jet quenching models is the 
initial color-charge density of the medium. By studying 
Cu+Cu collisions in the range y/SNN ~ 20 — 200 GeV 
this parameter can be varied with essentially no change 
in transverse size and shape of the reaction zone. More- 
over, the enhancement of hadron yields due to multiple 
soft scattering of the incoming partons ("nuclear Ax" or 
"Cronin enhancement") is expected to increase towards 
smaller y/s^N Mi thus the interplay between this en- 
hancement and the suppression due to parton energy-loss 
can be studied. 

In this letter we present invariant ir° yields for Cu+Cu 
collisions at y/sWW = 22.4,62.4, and 200 GeV. Refer- 
ence data for p+p collisions at yfs = 62.4 GeV and 
200 GeV were taken with the same experiment 13, 14j . 
At y/SNN = 22.4 GeV p+p reference data were obtained 
from a parameterization of the world's data on ir° pro- 
duction [lo| . 

Neutral pions were measured via their ir° — > 77 decay 
branch with the electromagnetic calorimeter (EMCal) of 
the PHENIX experiment [16[ . The EMCal comprises two 
calorimeter types: 6 sectors of a lead scintillator sam- 
pling calorimeter (PbSc) and 2 sectors of a lead glass 
Cherenkov calorimeter (PbGl). Each sector is located 
~ 5m from the beamline and subtends \rj\ < 0.35 in 
pseudorapidity and Aip — 22.5° in azimuth. Owing to 
the PbSc (PbGl) granularity of At] x Aip = 0.011 x 0.011 
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TABLE I: Data sets presented in this paper with the number 
of analyzed events. For the ERT triggered data the number 
of equivalent minimum bias events is given. 



system ^s NN 



St! 



TvrERT / 7vrsampled \ 
Jv evt l JV evt I 



Cu+Cu 22.4 GeV 75 - 90% 5.8 ■ 10 6 
Cu+Cu 62.4 GeV (88 ± 4)% 192 • 10 6 

Cu+Cu 200 GeV (94 ± 2)% 794 • 10 6 15.5 • 10 6 (4720 ■ 10 6 ) 



(0.008 x 0.008) the probability that the two photon show- 
ers from a ir° decay result in partially overlapping clus- 
ters is negligible up to a ir° pt of 12 GeV/c (15 GeV/c). 
The energy calibration of the EMCal was corroborated 
by the position of the 7r° invariant mass peak, the energy 
deposited by minimum ionizing charged particles travers- 
ing the EMCal (PbSc), and the correlation between the 
measured momenta of electron and positron tracks iden- 
tified by the ring-imaging Chercnkov detector and the 
associated energy deposited in the EMCal. These stud- 
ies showed that the accuracy of the energy measurement 
was better than 1.5%. 

The total number of analyzed Cu+Cu events for the 
three energies is shown in Table [T] The minimum bias 
(MB) trigger for all reaction systems was provided by 



< 



\V\ 



< 



Beam-Beam-Counters (BBC's) located at 3.0 
3.9. The reaction vertex along the beam axis, determined 
from the arrival time differences in the BBC's, was re- 
quired to be in the range \z\ < 30 cm. An additional 
hardware trigger (ERT) on high-px photons/electrons 
was employed in Cu+Cu at ^snn — 200 GeV. This trig- 
ger was based on the analog energy signal measured in 
overlapping 4x4 towers of the EMCal in coincidence with 
the MB trigger condition. The ERT reached a efficiency 
plateau for photon energies E > 4 GeV. 

The centrality selection in Cu+Cu at ^/snn = 
200 GeV and ^snn = 62 GeV was based on the 
charge signal of the BBC's which is proportional to the 
charged particle multiplicity in the respective pseudora- 
pidity range. The BBC trigger efficiency (etrig) for these 
systems was determined with the aid of the HIJING event 
generator and a full GEANT simulation of the BBC re- 
sponse (see Table H}. At ^snn = 22.4 GeV centrality 
classes were defined based on the charged particle mul- 
tiplicity measured with the pad chamber (PCI) detector 
(\rj\ < 0.35). The measured PCI multiplicity distribution 
was accurately reproduced in a Glauber Monte Carlo cal- 
culation [13] and centrality classes were determined by 
identical cuts on the measured and simulated PCI mul- 
tiplicities. The estimated BBC trigger efficiency given in 
Table [J results from a comparison of the simulated and 
the measured PCI multiplicity distributions. The results 
of the Glauber calculation 17| for Cu+Cu collisions at 
22.4, 62.4, and 200 GeV using inelastic cross sections of 
32.3, 35.6, and 42 mb, respectively, are given in Table HT1 



TABLE II: Glauber Monte Carlo calculations for Cu+Cu col- 
lisions at 22.4, 62.4, and 200 GeV. The N coll systematic un- 
certainty at 62.4 and 200 GeV is ~ 12%, almost independent 
of JVcoii. At 22.4 GeV the relative uncertainty of N co n can be 
parameterized as 0.094 + 0.173e" 0272JVc ° 11 . 





22.4 GeV 


62.4 GeV 


200 GeV 




(N A 




(N A 

\J ''part / 


(N n \ 




\J Vcoll/ 


0-10% 


92.2 


140.7 


93.3 


152.3 


98.2 


182.7 


10-20 % 


67.8 


93.3 


71.1 


105.5 


73.6 


121.1 


20-30 % 


48.3 


59.7 


51.3 


67.8 


53.0 


76.1 


30-40 % 


34.1 


38.0 


36.2 


42.6 


37.3 


47.1 


40-50 % 


23.1 


22.9 


24.9 


26.2 


25.4 


28.1 


50-60 % 


15.5 


13.9 


16.1 


15.0 


16.7 


16.2 


60-70 % 










10.4 


9.0 


70-80 % 










6.4 


1.9 


80-94 % 










3.6 


2.4 


60-88 % 






7.0 


5.5 







Neutral pions yields were measured on a statistical 
basis by calculating the invariant mass of all photon 
pairs in a given event and counting those within the 
7T° mass range. The background of combinatorial pairs 
was calculated by pairing photon hits from different 
events. Only photon pairs with an energy asymmetry 
\Ei - E 2 \/{E 1 + E 2 ) < 0.7 were accepted. The raw tt° 
yields were corrected for the geometrical acceptance and 
reconstruction efficiency. The latter takes into account 
the loss of 7r°'s due to photon identification cuts, the en- 
ergy asymmetry cut, inactive detector areas, and photon 
conversions. Moreover, it corrects the distortion of the 
7T° spectrum which results from the finite energy resolu- 
tion in conjunction with the steeply falling spectra and 
shower overlap effects. For Cu+Cu at ^snn = 200 GeV 
the transition between the minimum bias and the ERT 
sample occurs at p^ = 8 GeV/c. The final spectra were 
calculated as the weighted average of the PbSc and PbGl 
results, which agree well within the uncertainties. 

The main systematic uncertainties of the 7r° spec- 
tra result from the 7r° peak extraction, the reconstruc- 
tion efficiency, and the EMCal energy calibration. For 
Pt > 2 GeV/c the peak extraction uncertainty is ~ 4% 
for all systems, approximately independent of pt- The 
uncertainty in the reconstruction efficiency was estimated 
to be ~ 15 % for the three Cu+Cu analyses. It includes 
uncertainties of the photon identification cuts, the en- 
ergy resolution, and the modeling of shower overlap ef- 
fects. The uncertainty in the EMCal energy scale of 1.5 % 
translates into an uncertainty in the yields that increases 
from - 8 % at p T = 3 GeV/c to 15 % at p T = 6 GeV/c. 
The high-px part of the spectra in Cu+Cu at 200 GeV 
measured with the ERT trigger is subject to an additional 
uncertainty of 10% related to the ERT trigger efficiency 
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FIG. 1: Invariant it yields in central Cu+Cu collisions (a) and invariant 7r cross sections in p+p collisions (b) at ^/snn = 
22.4,62.4,200 GeV [II G3. The error bars represent the quadratic sum of the statistical and total systematic uncertainties. 
Plotted as a function of xt = 2px/v / s ( c ) the p+p data exhibit an approximate xt scaling. 



and normalization. 

PHENIX has not yet acquired a p+p data set at y/s = 
22.4 GeV. Therefore world data on charged and neutral 
pion production in the range 21.7 < ^/s < 23.8 GeV were 
scaled to y/s = 22.4 GeV and fit in the range < pt ^ 
7 GeV/c with Ed 3 a/dp 3 = A(e a ^ +h ) n (^ /2 - p T ) m 
where A = 1.22 • 1CT 17 mbGeV- 2 c 3 , a = 0.053 GeV^c, 
b = -0.884, n = -15.25, and m = 4.653 0. The scaling 
correction was determined with a next-to-leading-order 
QCD calculation. The scaling correction was largest for 
</s = 23.8 GeV and reduced these spectra by ~ 30% [l5j ]. 
The parameterization is consistent within ±25 % with the 
existing n° and ^ measurements without discernible pt- 
dependent systematic deviations. 

The 7T° pt spectra for p+p and central Cu+Cu colli- 
sions (0-10% of <tP^ +Cu ) at y/sWJi = 22.4, 62.4 [l|, and 
200 GeV @ are shown in Fig.EDi and[lj3. At sufficiently 
high pt where pion production in p+p collisions is dom- 
inated by fragmentation of jets, QCD predicts a scaling 
law v /j nrff ( XT 'V / s) Ed 3 a/dp 3 — G{xt) with a universal 
function G(xt) where xt = 2pt/\/~s |18|. Fig. [TJ; shows 
that such a scaling in xt is indeed observed for p+p colli- 
sions at 22.4, 62.4, and 200 GeV, consistent with previous 
observations [l9| . The xt values at which the universal 
curve G(xt) is reached indicate that particle production 
is dominated by hard processes for pt ^ 2 GeV/c for the 
three considered energies. 

Nuclear effects on high-px tt° production can be quan- 
tified with the nuclear modification factor 



-Raa(pt) 



(l/N c A f)d 2 N AA /d PT dy 
(Tab) x d 2 a pp /dpTdy 



(1) 



where (Tab) = (N co i\) / a p n p cl . In the absence of nuclear 
effects Raa = 1 for Pt 2 GeV/c where pions result 
from hard scattering processes. Raa (pt ) for the 0— 10% 
most central Cu+Cu collisions at 22.4, 62.4, and 200 GeV 



is shown in Fig. [2j The suppression at 62.4 GeV (Raa ~ 
0.6 for p T > 3 GeV/c) and 200 GeV (Raa ~ 0.5 - 0.6 
for pt ?i 3 GeV/c) is consistent with expectations from 
parton energy-loss. The Raa > 1 in Cu+Cu at 22.4 GeV 
is similar to the enhancement by a factor ^1.5 (at pt ~ 
3 GeV/c) observed in p+W relative top+Be collisions 
at y/sWlT = 19.4 GeV and 23.8 GeV ;2Q!]. For a similar 
number of participants the Raa in Cu+Cu at 22.4 GeV 
agrees with the Raa in Pb+Pb collisions at 17.3 GeV 

For pt ^ 3 GeV/c the measured nuclear modification 
factors at 62.4, and 200 GeV are consistent with a nu- 
merically evaluated parton energy-loss model described 
in |2l|, [22j as indicated by the comparison in Fig. [2] 



This calculation takes into account shadowing from co- 
herent final state interactions in nuclei [231, Cronin en- 



hancement [2J], initial state parton energy- loss in cold 
nuclear matter 25] , and final state parton energy-loss 



in dense partonic matter [9j, |2lJ, [22|. The Cronin en- 
hancement measured in p+A collisions is described well 
by this model [3]. The initial gluon rapidity density 
dN g /dy which characterizes the medium was not fit to 
the Raa values, but instead was constrained by mea- 
sured charged-particle multiplicities and the assumption 
of parton-hadron duality (dN 6 /dy = ndrj/dydN^/dr] 
with k = 3/2 ± 30% and drj/dy = 1.2 at all energies) 
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22j . The average fractional energy losses AE/E for 
a quark (gluon) with E — 6 GeV corresponding to the 
dNS/dy ranges in Fig. [U are 0.13 - 0.19 (0.29 - 0.42), 
0.16-0.20 (0.35-0.44), 0.20-0.28 (0.44-0.63) in central 
Cu+Cu collisions at 22.4, 62.4, and 200 GeV, respectively 
22]. For Cu+Cu at ^/s^a? = 22.4 GeV the calculation 
is also shown without final state parton energy-loss. The 
measurement is consistent with this calculation but does 
not rule out a scenario with parton energy-loss. 

Fig. [3] shows that the 7r° suppression in the range 
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-Vltev,\/s^;=22.4 GeV : 



1 1 1 1 1 1 r— 

22.4 GeV 
62.4 GeV 
200 GeV 

■ Vitev, 22.4 GeV, 130 < dN°/dy < 185 

□ Vilev, 62.4 GeV, 175 < dl^/dy < 255 

□ Vilev, 200 GeV, 255 < dN B /dy < 370 




10 



P T 1 (GeV/c) 







> 








\±> 




m 


1.5 


n 




V 








o!~ 




V 




m 


1 










A 




< 




< 




DC 




V 


0.5 




FIG. 2: Measured n° Raa as a function of pt for the — 10% 
most central Cu+Cu collisions at ^/sjvjv = 22.4, 62.4, 200 GeV 
in comparison to a jet quenching calculation (2ll. I22I]. The er- 
ror bars in this figure (and in Fig. [3} represent the quadratic 
sum of the statistical uncertainties and the point-to-point un- 
corrected and correlated systematic uncertainties. The boxes 
around unity indicate uncertainties related to (iV co ii) and ab- 
solute normalization. The bands for the theory calculation 
correspond to the assumed range of the initial gluon density 
dN s /dy. The thin solid line is a calculation without parton 
energy-loss for central Cu+Cu at -^/sjvjv = 22.4 GeV. 

2.5 < pt < 3.5 GeV/c increases towards more central 
Cu+Cu collisions for ^/JWW = 62.4, 200 GeV. On the 
other hand, Raa at y^s^F/v = 22.4 GeV remains approxi- 
mately constant as a function of N par t , suggesting either 
that the Cronin enhancement depends only weakly on 
centrality or that in this energy range parton energy-loss 
is offset by the larger effect of Cronin enhancement over a 
broad range of centrality. It appears from these data that 
in Cu+Cu collisions between t/snn = 22.4 and 62.4 GeV 
parton energy-loss will start to prevail over the Cronin 
enhancement, resulting in a net suppression. 

In summary, for the first time ir° spectra for the 
same nuclear colliding system (Cu+Cu) were measured 
in the same experimental setup over a wide range of ener- 
gies (^/sWn = 22.4,62.4, and 200 GeV). Nuclear effects 
were studied using measured p+p 7r° reference spectra 
from PHENIX at 62.4 and 200 GeV, and a parameteri- 
zation of world data at 22.4 GeV. High-px tt° yields in 
central Cu+Cu collisions at 62.4 GeV and 200 GeV are 
suppressed, suggesting that parton energy-loss is a signif- 
icant effect in these systems. At 22.4 GeV tt yields for 
Pt > 2 GeV/c are not suppressed. The Raa measured in 
central Cu+Cu at 22.4 GeV is consistent with Cronin en- 
hancement alone but does not rule out parton energy-loss 
effects. The measurements of high-px tt production over 
a factor ~ 10 in center-of-mass energy presented in this 
letter provide a unique constraint for jet-quenching mod- 
els and demonstrate that parton energy-loss starts to pre- 



FIG. 3: The average Raa in the interval 2.5 < pt < 
3.5 GeV/c as a function of centrality for Cu+Cu collisions 
at ^/JWn = 22.4, 62.4, and 200 GeV. The shaded bands rep- 
resent jet quenching calculations at three discrete centralities 
(AW ~ 10, 50, 100) [U, HI- The boxes around unity repre- 
sent the normalization and (iVcoll) uncertainties for a typical 
iVcoii uncertainty of 12 %. 



vail over the Cronin enhancement between ^/snn = 22.4 
and 62.4 GeV. 
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